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Abstract 



We prove the integrability of the short pulse equation derived recently by Schafer and Wayne 
from a hamiltonian point of view. We give its bi-hamiltonian structure and show how the recursion 
operator defined by the hamiltonian operators is connected with the one obtained by Sakovich and 
Sakovich. An alternative zero-curvature formulation is also given. 
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The cubic nonlinear Schrodinger (NLS) equation is used to describe the propagation of optical 
pulses in nonlinear media. Extending the work of Altermnan and Rauch ^ the short pulse (SP) 
equation 

%ct = U + - {U 3 ) xx , (1) 

was derived by Schafer and Wayne as an alternative equation to the NLS equation to describe 
the evolution of very short optical pulses in nonlinear media. Also, they made numeric studies that 
proved that as the pulse length gets short the SP equation is a much better approximation to the 
solution of Maxwell's equation than the NLS equation These results support the integrability 
of the SP equation based on numerical techniques. In an interesting paper Sakovich and Sakovich 
[3] studied the integrability of the SP equation from a Lax or zero-curvature point of view. They 
found that the equation Q possesses a Lax pair of the Wadati-Konno-Ichikawa type 15 and that 
through a chain of transformations this equation can be related to the sine-Gordon equation. This 
result sets the equation Q as ultrashort pulse integrable alternative to the NLS equation. Also, 
using cyclic basis techniques these authors obtained the following recursion operator 

R = d 2 —d{l + ul)^d- \ y (2) 

y-xx (i + uiyi 1 

It is well known that the integrability of a nonlinear equation can be approached from a Lax 
representation. Using a spectral transform method analytic solutions can in principle be obtained. 
However, these calculations turn out to be very complex and an alternative approach to integrability 
is desirable. An algebraic approach based on a hamiltonian formulation for infinite dimensional 
dynamical systems is very powerful and was intensely developed over the last several years after 
the seminal works of Gardner Zakharov and Faddeev [7] and Magri [SJ. In this letter we will 
show that the SP equation is bi-hamiltonian. These two compatible hamiltonian operators allow 
us to establish the integrability of the SP equation from a hamiltonian point of view without the 
introduction of Backlund transformations or the use of Lax equations. As it is well known [S] 
from these structures we can construct a recursion operator that yields a hierarchy of conserved 
equations and charges. For the SP equation its bi-hamiltonian structure will result in a recursion 
operator whose inverse is exactly (|2"|). A recursion operator (bi-hamiltonian structure) can also 
provide insights on a Lax or zero curvature formulation for an integrable system. Therefore, we 
will also obtain an alternative zero-curvature formulation for the SP equation which will provide 
in a natural way the hierarchy of charges and equations associated with it. 

The equation (^Q) studied in Ref. [I] is related to the original form used in [21 E] by scale 
transformations of the dependent and independent variables, but in order to obtain a hamiltonian 
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description for the SP equation let us integrate (^Q) once with respect to x to obtain the nonlocal 
equation 

u t = (<9 _1 u) + - u 2 u x . (3) 

Nonlocal equations, as well hierarchies of equations, of this sort have been studied extensively in 
the last years 

Introducing the Clebsch potential 

u = <p x , (4) 

the equation Q can be obtained from a variational principle, 5 J dtdx C, with the first order 
lagrangian density 

c = \ Mx ~ hf* + \^ ■ (5) 

Using the Dirac's theory of constraints ^U] we can obtain the hamiltonian and the hamiltonian 
operator associated with © 



v x = d, 



H 2 = I dx 



4 l/o-l \2 

—u (o u) 

24 2^ ; 



(6) 



It can be easily checked that H 2 is conserved by the SP equation. In this way the SP equation in 
the form (j3J) can be written in hamiltonian form as 

ou 

As it is well known a system can be shown to be integrable if it is bi-hamiltonian [HI II 1| I12j , 
i.e., the system has a hamiltonian description with two distinct hamiltonian structures that are 
compatible. In fact, as can be easily checked, the SP equation can also be written as 

OU 

with 

T> 2 = d~ l + u x d~ l u x , 

H^^fdxu 2 . (7) 



3 



Again, Hi is conserved by the SP equation. Also, the skew symmetry of this hamiltonian structure 
is manifest. The proof of Jacobi identity for this structure as well the compatibility with the 
hamiltonian structure in ((HJ) can be determined through the method of prolongation 11 in a 
straightforward way. Therefore, the bi-hamiltonian structures © and (J7J) also guarantees the 
integrability of the SP equation from a hamiltonian point of view and also provide us with a 
natural recursion operator defined by 

R = V 2 V^ (8) 

or 

R= (l + u 2 x -u x d- l u xx )d- 2 . (9) 
It can immediately be checked that the inverse of R is exactly R obtained by Sakovich and Sacovich 

From the recursion operator (JSJ) a whole hierarchy of equations and conserved charges can be 
constructed and the details will be published elsewhere JH]- However, since the system is bi- 
hamiltonian the SP equations and the other equations of the hierarchy will have the general form 

u t = (V 2 X) - X (DxX) , (10) 

where A represents a spectral parameter and X is an arbitrary function of the dynamical variable 
as well as of the spectral parameter. This equation can alternatively be obtained from the zero 
curvature condition (see ^1] and references therein) 

d t Ai - d x A - [A , Ai] = . (11) 

From the SP equation the canonical dimensions of the variables are given by [u] = —1, [x] = — 1, 
and [t] = 1. Also, from (J10[) we have [X] = —1 and [A] = —2. Let us choose, taking into account 
(|ll(l . the gauge fields with matrices 

] = -l [] = 1 

A = I I , Ai = I 1 (12) 




] = -3 [] = -! 

where [ ] represents the dimensionality of the matrix element. Choosing the gauge fields as 

-A- 1 /2(a-i u ) A -V2 

(d- 2 x) + u{d- l u x x) - {d- l uu x x) - x-^id^u) 2 - xx a- 1 / 2 ^- 1 ^) 



A 



o 
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A" 1 / 2 

Ai = | | , (13) 

u A-V2 

the zero curvature condition (|TT|) indeed yields the dynamical equations in ifTUjl 

ut= (d~ l +u x d- l u x )X-\dX . (14) 

Dynamical variables are independent of the spectral parameter, therefore, the function X must 
depend on A for this equation to be meaningful. Making a Taylor expansion in A of the form 

n 

X(X,u) = J2^~ j Xj(u), (15) 

3=0 

with 

x ° = TfTT Ws ' (16) 

and substituting (fT3)) and (fTBT) into ((Til) , we obtain 

V 2 X j {u) = V x X j+x {u) , j = 0, 1, 2, . . . , n - 1 , 

n tn = P 2 X n . (17) 

If we identify 

then (JUJ) gives the dynamical equations of the hierarchy (for any n) , the two hamiltonian structures 
of the system as well as the recursion relations between the conserved charges of the hierarchy. For 
completeness, the first flows are 



u to = u x , 

u tl = (<9 _1 -u) + \u 2 u x , 



(19) 



u t2 = (d- 3 u) + i(<rV) + u x (d- 1 (u x {d- 2 u))) + 1 u\ x , 



The n = flow is the chiral equation and the n = 1 flow is just the SP equation (j3J). The conserved 
charges for the hierarchy can, of course, be determined recursively from (|17|) and give the infinite 
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set of conserved hamiltonians 



Hq = — J dx y / l _ +~uf , 



dxu 2 , 



H 2 
H 3 



dx 



dx 



1 4 I/q-I \2 
— -11 [O U) 

24 2 V ' 



— u +-(9 u) +-(d u)u--(d u) u 



(20) 



Our zero-curvature formulation (|13|) for X = u and A = is nonlocal and yields the SP equation 
in the form (jSJ). Sakovich and Sakovich have also obtained a zero-curvature for the SP equation 
written in the local form (fl|). Both zero-curvatures are related by some sort of gauge transformation. 

In this letter we have supported the integrability of the short pulse equation through its bi- 
hamiltonian structure. We related our work with the previous work of Sakovich and Sakovich [1] 
comparing their recursion operator with ours. A zero-curvature formulation providing a hierarchy 
of conserved charges and equations, containing the SP equation, was given. Following the results 
of [2] this hierarchy can be extended to include other equations such as the elastic beams under 
tension equation ^3] and the details can be found in |13j . 

Acknowledgments 

This work was supported by CNPq (Brazil). 

References 

[1] D. Alterman and J. Rauch, Phys. Lett. A 264, 390 (2000). 

[2] T. Schafer and C. E. Wayne, Physica D 196, 90 (2004). 

[3] Y. Chung, C. K. R. T. Jones, T. Shafer and C. E. Wayne, Nonlinearity 18, 1351 (2005). 

[4] A. Sakovich and S. Sakovich, J. Phys. Soc. Jpn. 74, 239 (2005). 



6 



[5] M. Wadati, K. Konno and Y. H. Ichikawa, J. Phys. Soc. Jpn. 47, 1698 (1979). 
[6] C. S. Gardner, J. Math. Phys. 12, 1548 (1971). 

[7] V. E. Zakharov and L. D. Faddeev, Func. Anal. Appl. 5, 280 (1971). 
[8] F. Magri, J. Math. Phys. 19, 1156 (1978). 

[9] J. C. Brunelli and G. A. T. F. da Costa, J. Math. Phys. 43, 6116 (2002); J. C. Brunelli, A. 
Das and Z. Popowicz, J. Math. Phys. 45, 2646 (2004); J. C. Brunelli and A. Das, J. Math. 
Phys. 45, 2633 (2004). 

[10] P. A. M. Dirac, Lectures on Quantum Mechanics, Belfer Graduate School of Science Mono- 
graphs, vol. 2 (New York, 1964); K. Sundermeyer, Constrained Dynamics, Lecture Notes in 
Physics, vol. 169 (Springer, Berlin, 1982). 

[11] P. J. Olver, Applications of Lie Groups to Differential Equations, 2nd ed. (Springer, Berlin, 
1993). 

[12] M. Blaszak, Multi-Hamiltonian Theory of Dynamical Systems (Springer, Berlin, 1998). 

[13] J. C. Brunelli, "The Short Pulse Hierarchy", to appear in the Journal of Mathematical Physics. 

[14] A. Das and S. Roy, Mod. Phys. Lett. All, 1317 (1996). 

[15] Y. Ichikawa, K. Konno and M. Wadati, J. Phys. Soc. Jpn. 50, 1799 (1981). 



7 



